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Abstract. This paper presents the evolution of the clustering of the main population of galaxies from z ~ 2 to 
2 = 0.2, from the first epoch VIMOS VLT Deep Survey (VVDS), a magnitude limited sample with 17.5 < Iab < 
24. The sample allows a direct estimate of evolution from within the same survey over the time base sampled. We 
have computed the correlation functions ^(rp,7r) and Wp{rp), and the correlation length ro{z), for the VVDS-02h 
and VVDS-CDFS fields, for a total of 7155 galaxies in a 0.61 deg^ area. We find that the correlation length in 
this sample slightly increases from z — 0.5 to z = 1.1, with ro{z) = 2.2 — 2.9 h~^ Mpc (comoving), for galaxies 
comparable in luminosity to the local 2dFGRS and SDSS samples, indicating that the amplitude of the correlation 
function was ~ 2.5 times lower at z ~ 1 than observed locally. The correlation length in our lowest redshift bin 
z = [0.2,0.5] is ro — 2.2 h~^ Mpc, lower than for any other population at the same redshift, indicating the low 
clustering of very low luminosity galaxies, 1.5 magnitudes fainter than in the 2dFGRS or SDSS. The correlation 
length increases to ro ~ 3.6 h~^ Mpc at higher redshifts z = [1.3,2.1], as we are observing increasingly brighter 
galaxies, comparable to galaxies with M b a a = —20.5 lo cally. We compare our measurement to the DEEP2 



measurements in the range z = [0.7, 1.35] I Coil et all 12004^ and find comparable results when applying the same 



magnitude and color selection criteria as in their survey. The slowly varying clustering of VVDS galaxies as 
redshift increases is markedly different from the predicted evolution of the clustering of dark matter, indicating 
that bright galaxies traced higher density peaks when the large scale structures were emerging from the dark 
matter distribution 9 — 10 billion years ago, being supporting evidence for a strong evolution of the galaxy vs. 
dark matter bias. 



Key words. Cosmology: observations - Cosmology: deep redshift surveys - Galaxies: evolution - Cosmology: large 
scale structure of universe 



1. Introduction 

The evolution of the clustering of galaxies is a key diag- 
nostic element to test the evolution of the universe, al- 
lowing direct comparison between observations and the- 
ory. In the current paradigm of galaxy formation and evo- 
lution, dark matter halos that contain galaxies are ex- 
pected to merge and grow under the action of gravity. 
This translates into a continuous evolution of the correla- 
tion function ^(r, z) of dark matter halos, now quite well 
understood fr om extensive high reso l ution numerical sim- 
ulations (see IWeinberg et all '20 oJ iBenson et all l200lt 
ISomerville et all l200lt Ixauffman n et all ll999|r As a di- 
rect measurement of the space distribution of dark matter 
halos is not yet feasible, we are compelled to use galaxies 
as indirect tracers of the dark matter. Unfortunately, as 
galaxies are complex physical systems, their relationship 
to the underlying mass, the "bias", is difficult to estimate. 
As galaxies and dark matter evolve, the bias may evolve 
and relating the measurements of galaxy clustering to the 
evolution of the total mass is not easy, with the bias shown 
to dep end upon galaxy type , luminosity and local environ- 
ment l)Norberg et all l2002|) . 

The most straigthforward indicator of galaxy cluster- 
ing is the correlation function ^(r), representing the excess 
probability over random of finding a galaxy in a given vol- 
ume, at a fixed distance from another galaxy. The shape 
and amplitude of the galaxy correlation function at the 
current epoch is now established to high accuracy, ^(r) is 
well described by a power law f (r) = (r/rn)"''' over scales 
0.1 - inh-i Mnc llDavis Peeble4 Fms3 iHawkins et all 
l200.'4 IZehavi et all EoM) . with a more refined modelling 
requiring some extra power over thi s shape for separa- 
tions larger than 2 - 3 h~^ Mpc ijGuzzo et all Il99lt 
IZehavi et al'll2004j) . a feature possibly encoding informa- 
tion on the relation between galaxies and their host dark- 
matter halos. The local clustering measurements have 
shown that the correlation length tq increases from late 
types to early type galaxies, from low luminosity to high 
luminosity a nd from low to high galaxy density envi- 
ronments (se e Giova nelli et allll98(iHBenoist et allll99(it 
IGuzzo et all .1997.: iNorberg et all l200ll). with lum inos- 
ity being the dominant effect l|Norberg et al 1 120021) . The 
most recent estimates of the correlation length from the 
2dFGRS and SDSS vary from rg = 3 for late type star- 
forming galaxies in low density environments to tq = 5 — 6 
for galaxies with Af* = —19.5, with the clustering ampli- 
tude reaching Tq = 7.5 h~^ Mpc for galaxies four times 

At higher redshifts the situation is less clear. Analysis 
of the projecte d angular correl a tion function fav ors a sta- 
ble clustering l|Postman et al 1 119981: iRoche et al.. 1999 : 
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iMcCracken et all l200lt ICabanac et"ail l2000() but it re- 
quires a priori knowledge of the redshift distribution of 
the galaxy population sampled. A variety of results have 
been obtained from smaller spectroscopic samples, with 
comoving correl ation lengths rn in the range 2 — 5 h~^ 
Mpc at z ~ 0.5 jiLe Fevre et al.Ul99 6: Small et~all Il999t 
IShepherd et"ani200lHCarlberg et al.L.1999} andrn = 3-5 
h~MVIpc at z ^ 3 ijGiavalisco et allll998l:lFoucaud et all 
Recently, first results from the DEEP 2 survey have 
been presented, indicating a correlation length tq = 3.53± 
0.81 h-i Mpc in z=[0. 7,0.91. and r n = 3.12±0.72 h^^ Mpc 
in z=[0.9,1.35] l)Coil et all l2004|) . The main difficulty in 
interpreting these results in terms of evolution of the clus- 
tering is to relate the population of galaxies observed at 
a given redshift to a lower redshift population of "descen- 
dants" , identified from a well-defined selection function 
enabling comparisons. Some of these surveys are targeted 
to specific classes of galaxies, pre-selected via photomet- 
ric methods, whose relation to the global population is 
not obvious. The most notable example is represented by 
galaxies selected vi a the Lyman-break technique around 
z ~ 3 and z ^ 4 ijSteidel et all Il998l) . which display a 
clustering strength similar to present-day normal galaxies 
and therefore represent a very biased po pulation, possibly 
the p recursors of giant cluster ellipticals l|Governato et all 

On the other hand, even when selecting purely 
magnitude-limited samples, one cannot avoid being af- 
fected by the complex dependence of clustering on mor- 
phology and luminosity evidenced by the wide range of 
correlation lengths measured in the local Universe, and 
its evolution as a function of redshift. One needs to ob- 
serve samples at increasingly high redshifts with the same 
luminosities, colors (type) and local environments in or- 
der to derive the evolution of the clustering of galaxies and 
hence attempt to derive how the correlation properties of 
the mass evolve. At high redshifts, the natural observa- 
tional bias is to sample increasingly brighter and more 
actively star forming galaxies, which may have a direct 
impact on our current vision of the evolution. 

Finally, high-redshift samples of spectroscopically 
measured galaxies have been inevitably limited so far to 
relatively bright objects in small areas on the sky, which 
contributes to increase the scatter between independent 
measurements, further complicating their interpretation. 

In this paper we present the first attempt to mea- 
sure the evolution of the clustering in a consistant way 
accross the redshift range 0.2 < z < 2.1, using 7155 
galaxies from the VIMOS VLT Deep Survey (VVDS) over 
more than 0.61 deg^. The VVDS is designed to sample 
the high redshift population of galaxies in the most unbi- 
ased way possible, using a simple magnitude selection in 
the range 17.5 < Iab < 24, using severa l indep endant 
fields up to 4 deg2 each dLe Fevre et all l2005a|) . This 
analysis uses the high-quality First Epoch WDS sam- 
ple, which includes 6117 galaxy redshifts in the VVDS-02h 

T — — 1 ' y 

IjLe Fevre et al., 2005a) and another 1368 galaxies in the 
WDS-CDFS (,Le Fevre et a,1...2004bil field. We have mea- 
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sured the redshift-space correlation function f (rp, tt), and 
computed the projected function Wp{rp), to recover the 
value of the galaxy correlation length ro{z) up to z ~ 2, 
therefore tracing the evolution of the clustering over more 
than 10 Gy, or 70% of the current age of the universe. 

In section 2 we recall the properties of the VVDS First 
Epoch sample. In section 3, we describe how the correla- 
tion function has been comput ed, referring in lar ge part 
to the accompanying paper by IPollo et al.l l(200,'i|) which 
describes all the methods set up to validate the measure- 
ments and compute the errors. In section 4, we present the 
results in terms of the evolution of tqIz), and we compare 
our results to previous measurements where possible. In 
section 5, we discuss the evolution of the clustering of the 
global population of galaxies from z ~ 2, before conclud- 
ing in section 6. 

This paper is the first in a series of papers to study the 
clustering of galaxies at high redshift from the VVDS first 
epoch data. Guzzo et al. (2005, in prep.) will present the 
clustering evolution from volume limited samples and in- 
fer the dependence of clustering upon luminosity; Meneux 
et al. (2005, in prep.) will present the differences in clus- 
tering observed as a function of galaxy type and its evo- 
lution and Polio et al. (2005, in prep.) will investigate the 
dependence and ev olution of cl ustering as a function of the 
local environment. iMarinoni et al. (200^ and Le Fevre et 
al. (2005, in prep.) will look at the evolution of the galaxy 
- dark matter bias, and subsequent papers will study the 
clustering from the redshift population 2 < z < 5. 

We have used a Concordance Cosmology with fi™ = 
0.3, and SIa — 0.7 throughout this paper. The Hubble 
constant is normally parameterized via h=Ho/100, to ease 
comparison to previous works, while a value Ho = 70 km 
s~^ Mpc~^ has been used when computing absolute mag- 
nitudes. All correlation length values are quoted in comov- 
ing coordinates. 

2. VVDS first epoch data 

2.1. The sample 

The VVDS-Deep sample is stricly selected in magni- 
tude in the range 17.5 < Iar < 24, from a com- 
plete deep photometri c survey l)Le Fevre et al.l l2004at 
iMcCracken et al ll200 ^ without any color or shape selec- 
tion l|Le Fevre et al.L l2004b'l. We have analysed two fields, 
the VVDS-02h and the VVDS-CDFS. Over the ~ 9600 
redshifts measured in these two fields, we have used in 
the following analysis only objects with a redshift confi- 
dence level higher th an or equal to c ± 80% ( quality flags 
2 to 9 as defined in iLe Fevre et alJ l 2004b() \ excluding 
QSOs. We will only mention briefly below the effect on 
measured correlations of relaxing the quality threshold, 
including the poorest redshift measurements (flag 1). The 
complete sample analysed concerns a total of 7155 galaxies 
in 2203 arcmin^, with 6137 galaxies in the 1750 arcmin^ 
VVDS-02h field with 0.2 < z < 2.1, and 1038 galaxies 
in the 453 arcmin^ VVDS-CDFS (Chandra Deep Field 
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Fig. 1. The redshift distribution in the VVDS-02h field. 
The filled histogram contains 6137 galaxies with quality 
flag > 2, and another 1093 galaxies with quality flag — 1 
(open dashed histogram), in the redshift range z=[0.2,2.1] 



South, ijGiacconi et all kOO^ ) area with 0.4 < z < 1.5. 
The accuracy of the red shift measurements is ~ 275 km/s 
l|Le Fevre et all l2004a|) . The distribution of redshifts for 
the VVDS-02h field is presented in Figure H The redshift 
distribution peaks at z ~ 0.8, and there are 300 galaxies 
with 1.3 < z < 1.5, and 132 galaxies with 1.5 < z < 2.1. 

2.2. The galaxy population mix 

We have split the sample into 6 redshift bins, as described 
in Table ni The rest frame B — Iab{0) color and absolute 
magnitude Mbab distribution within each bin is shown in 
Figure 121 from z = 0.2 to z = 2.1, and the mean values 
are reported in Tabled B — Iab{^) and Mb^b have been 
computed using template fitting of the photometric spec- 
tral energy distribution in B,V,R, and I ba nds, to derive 
the fc(z) corrections (see lllbert et al.L [20051 for details). 

Up to redshift z = 1.3, the B — /ab(0) color distri- 
bution stays quite similar with increasing redshift, from 
blue star forming B — /ab(0) = to red B — Iab{0) ~ 3, 
while for 1.3 < z < 2.1, the reddest B - Iab{0) > 2.5 
galaxies, if present, are not observed. The magnitude se- 
lection of the VVDS-Deep, therefore, allows us to sample 
the global population of galaxies, for all galaxy types from 
late to early types, up to z ~ 1.3, while for z > 1.3, as 
we are selecting galaxies from their UV rest frame contin- 
uum at these redshifts, the VVDS is increasingly biased 
towards late-type, higher star-forming galaxies. We there- 
fore expect that in the farthest redshift bin, the clustering 
measurement in the VVDS is the result of the effects of 
looking at intrinsicaly more luminous and more actively 
star forming galaxies. 
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The range of absolute Mb^b magnitudes sampled is 
quite large, and changes strongly with redshift as shown 
in Figure 121 While at redshift z ~ 0.5, the absolute mag- 
nitude range sampled is —22 < Mb^b ^ ~16, only bright 
galaxies with Mbab — ~20 are sampled ai z > 1.3. The 
absolute luminosity of galaxies in the VVDS is shown 
to increase as redshif t increases, a,s com puted from the 
Luminosity Function l|llbert et alil2005h . If one assumes 
a pure luminosity evolution, galaxies observed at high red- 
shift are expected to have faded to fainter absolute mag- 
nitudes at the present time, by as much as Mb = 1.5 to 
2 magnitudes. This should be taken into account when 
comparing the clustering of the high redshift population 
to local populations. 

The consequences on the correlation function measure- 
ments of the change in the population sampled as a func- 
tion of redshift are discussed in Section 5. 



3. Computing the real-space correlation and 
correlation length ro(^) 

The methods apphed on the VVDS first epoch data to 
derive the real space correlation parameters a re described 
extens ively in the accompanying paper by IPollo et alJ 
(HqQ^. We summarize below the main elements of this 
method. 

3.1. Estimating correlation functions from the VVDS 

To measure the galaxy real-space correlation length Tq and 
slope 7 from our survey, we have used the projection of the 
bi-dimensional correlation function ffr^ ,7r) . This func- 
tion w as estimated using the well-known Landv & Szalavl 
119931) estimator 



No RRirp,7T) 



where Nq is the mean galaxy density (or, equivalently, 
the total number of objects) in the survey; Nji is the mean 
density of a catalogue of random points distributed within 
the same volume of the considered redshift bin; GG{r) is 
the number of independent galaxy-galaxy pairs with sepa- 
ration between r and r + dr; RR{r) is the number of inde- 
pendent random-random pairs and GR{r) the number of 
cross galaxy-random pairs within the same interval of sep- 
arations . A total of ^ 40, 000 random points have been 
used in each redshift bin, guaranteeing a sufficient den- 
sity to avoid shot-noise effects on small scales. The ran- 
dom sample follows exactly the same geometry, redshift 
distribution and observational pattern as in the galaxy 
data, while a specific weighting scheme is used to overcome 
the biases introduced by the slit-positioning software and 
other selection effects. These techniques have been exten- 
sively tested on a large number of mock VVDS surveys 
and have been shown to be able to recover the correct 
Wp{rp) correlation function to better than 10%, reducing 
any systematic effect to less than 5%. We have not applied 
this small correction to the data presented in this paper. 



This is discussed in detail in the accompanying paper by 
IPollo et alJ ll2005l) . 

Since we are not computing the correlation function 
from the full magnitude-limited survey altogether, there is 
no point here in using the so-called J3 minimum-variance 
weighting. This usually has been ad opted in the analy- 
sis of large fl ux-limited local surveys l|Fisher et all Il994t 
IGuzzo et al.l lioOO'). in which the sampling of the cluster- 
ing process varies dramatically between the nearby and 
distant parts of the sample. Its main scope is to avoid ex- 
cessive weighting of the most distant parts of the sample, 
where only sparse bright galaxies trace structures. Within 
each of our redshift bins, the density of objects varies only 
slightly and equ al weighting of the pairs is the most ap- 
propriate choice ( Fisher et all Il994l) . 



3.2. ^(rp,7r), Wp{rp) and the correlation length tq 

We have computed the two point correlation function 
^(rp,7r) in increasing redshift bins, selecting the bin 
boundaries to maximize the number of objects, hence the 
signal to noise of the correlation function, in each of the 
bins. 

The formalism developed by iDavis fc Peeblesl l)l983l) 
has been used to derive the real space correlation func- 
tion in the presence of redshift-space distortions along the 
line of sight. We integrate C(rp,7r) along the line of sight, 
computing 



Wp{rp) 



dy£.[rp,TT) 



irl + y^^' ,(2) 



where in practice the upper integration limit has to be 
chosen finite, to avoid adding noise to the result. After a 
few experiments, we used a value of 20 h~^ Mpc. 
If we assume a power-law form for ^(r), i.e. 



CM = - 



Wp{rp) can be written as 



Wpirp) 



-r(i)r(V) 



(3) 



(4) 



where F is the Euler Gamma function. 

Fitting the 'Wp{rp) measurements in each redshift bin 
then provides a measurement of ro(2:) and 7(2;). 

3.3. Computing errors 

The uncertainty associated with the computation of 
rQ{z) and 7(2) is largely dominated by cosmic variance. 
Although our spectroscopic sample is the largest avail- 
able to date at the redshifts probed, both in the num- 
ber of galaxies and area surveyed, only two fields have 
been sampled, with the VVDS-02h largely dominating 
over the VVDS-CDFS in terms of the number of galax- 
ies observed and area covered, and it is therefore inap- 
propriate to estimate errors on Wp(rp), ro(z), and 7(2;) 
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Fig. 2. The rest frame B — Tab color (left) and absolute magnitude Mb^b (right) from the VVDS-02h data up to 
z - 2. 



directly from field-to-field variations within our data set. 
Instead, we have been using ensemble errors derived from 
the scatter in the correlation function computed from 50 
mock VVD S surveys constru cted using the GallCS sim- 
ulations (Bl aizot et al.U2004l). as thoroug hly described in 
the parallel paper by l)Pollo et alJ. bOOsj) . These realistic 
mock samples specifically include the same selection func- 
tion and observational biases present in the VVDS First 
Epoch data set, and allow us to compute ensemble statis- 
tical averages and scatters of the quantities we measure 
on the real data. In particular, the standard deviation of 
the measured Wp{rp) from the mock VVDS samples pro- 
vides us with a realistic set of error bars that include also 
the effect of field-to-field variations due to fluctuations on 
scales larger than the observed field, i.e. cosmic variance 
(obviously assuming that the A-CDM power spectrum and 
normalization are a good match to the observed one, which 
is indeed the case). 



The comparison of the VVDS-02h and VVDS-CDFS 
correlation functions gives us an external check to this 
procedure, with a relatively noisy indication of the ampli- 
tude of cosmic variance, as described below. 

4. Results 

4.1. The VVDS correlation function from a magnitude 
limited sample 17.5 < Iab < 24 out to z ~ 2 

We have computed the correlation function ^{rp,n) and 
its projection Wp{rp) on both the VVDS-02h and VVDS- 
CDFS fields. 

The correlation function Wp{rp) is presented in the left 
panel of Figure 01 for the VVDS-02h field, and in Figure 
El for the VVDS-CDFS field. The number of galaxies ob- 
served in each field is indicated for each redshift bin in 
Figure so far the largest sample of galaxies used to 
compute the correlation function at these redshifts. The 
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large sample allows to compute Wp{rp) in 6 rcdshift bins 
up to z = 2.1 for the VVDS-02h field, and in 3 bins up 
to z = 1.5 for the VVDS-CDFS field. Error bars are en- 
semble errors computed as described in § 13.31 A positive 
correlation signal is measured out to at least 30 h^^ Mpc 
in all bins, and Wp{rp) is well described by a power law in 
the range 0.1 < rp < 10 h~^ Mpc (note however that any 
redshift space feature in ^(r) is smoothed out in Wp{rp) , 
which is its integral). The measured correlation function 
amplitude is relatively low at low redshifts z < 0.5, and 
stays essentially constant as a function of redshift, with 
just a mild increase in the farthest bins. A possible inter- 
pretation of these results is discussed in the next sections. 



4.2. Evolution of the galaxy real-space correlation 
function 

The measured values of the correlation length rQ{z) and 
the correlation function slope 7(2) computed from the fit- 
ting of Wp{rp) are reported in Tabled For the VVDS-02h 
field, we have used all Wp points for 0.1 < rp < 10 h^^ 
Mpc. For the VVDS-CDFS field, we have used points in 
the range 0.1 < Vp < 3 Mpc only, because of the 
smaller field size. We also report in Tabled the values of 
ro{z) obtained after fitting Wp with the slope 7 fixed to 
the average slope measured in the range z=[0.2,1.3]. In the 
VVDS-CDFS, a strong wall- like large scale structure has 
been identified at z ^ 0.735, wit h more than 130 galaxie s 
in a velocity range ±2000 km/s ()Le Fevre et al.U2004bj) . 
and is expected to strongly affect the correlation func- 
tion computation. The correlation function in this field has 
been computed in the redshift interval z= [0.6,0.8] which 
includes this strong structure, as well as in z=[0.4,0.7] 
and z=[0.8,1.5] to specifically exclude it, as reported in 
Tabled While we find a slightly higher correlation length 
in the VVDS-CDFS, the measurements are compatible, 
within the errors, with the values reported for the VVDS- 
02h. 

We have investigated how our results change when us- 
ing the full set of available redshifts, i.e. including also 
the 1300 less accurate meas urements (flag 1). whi ch are 
shown to be ~55% correct l)Le Fevre et al.L l2005ah . The 
result is that the measured correlation lengths are low- 
ered by ~ 3 — 5%, which we interpret as a consequence of 
the significant fraction of poorly measured redshifts which 
dilute the actual correlation function. 

The results strikingly show a slowly increasing corre- 
lation length over the complete redshift range z= [0.2,2.1]. 
The lowest value measured is in the lowest redshift bin 
probed, then tq is constant tq ~ 2.8 h~^ Mpc over the 
range z=[0.5,l.l] and increases slightly to tq ~ 3.6 h~^ 
Mpc over z=[l. 1,2.1]. When we fit the slope 7 of the cor- 
relation function at the same time as ro, it varies between 
1.67 and 1.96, slightly increasing with redshift, as reported 
in Tabled The average slope in the range z=[0.2,1.3] is 
7 = 1.76. 



4.3. Comparison with other surveys 

In redshift surveys of the local Universe (z < 0.2), the 
lowest values for the correlation length have been mea- 
sured for late-type star form ing galaxies as e.g . H-II 
galaxies (ro = 2.7 h"^ Mpc, llovino et aD ^M)) and 
infrared-sele c ted I RAS galaxies (ro = 3.76 h~^ Mpc, 
iFisher et al ](|l99l). Here we measure an even lower cor- 
relation length for galaxies with z=[0.2,0.5] in the VVDS, 
coherent with their very low intrinsic luminosity as dis- 
cussed in Sectional 

In order to compare our measurements at higher 
redshifts to the measurements of the DEEP2 survey 
l|Coil et al.Ll2004() . we have restricted our data by applying 
a-posteriori the same color and magnitude selection func- 
tion as the DEEP2 survey applied a-priori to pre-select 
its spectroscopic targets in the redshift range 0.7 — 1.35. 
With this "DEEP2" selection function {Rab < 24.1, 
B-R< 2.35(i? - /) - 0.45, R-I > 1.15 or B - R < 0.5) 
only 59% of the VVDS Tband magnitude-limited sample 
galaxies are selected (in the desired redshift range). The 
correlation function Wp(rp), and the corresponding ro(z) 
and 7(2) of this sample are shown in FigureEl We find that 
ro = 3.51 ±0.63 h'^ Mpc for the full z=[0 7,1.35] sample, 
compared to 3.19±0.51 bv lCoil et al. I 1I2OO4I) . Separating in 
the same redshift bins, and setting the slope to 7 = 1.66 
as measured in DEEP2, we find that ro = 3.37 ± 0.67 
h-i Mpc in z^ [0. 7,0.9], rp = 3.05 ± 0.47 h'^ Mpc in 
z=[0.9,1.35], while ICoil etaP ll2004l) find ro = 3.53 ± 0.81 
h-i Mpc and ro = 3.12 ± 0.72 h'^ Mpc, respectively, as 
shown in FigureEl Given the relative uncertainties of both 
surveys, these results are therefore in excellent agreement. 
The correlation function for the DEEP2 survey therefore 
differs from the correlation function of the complete WDS 
sample because of the different selection function applied 
that selects galaxies brighter than in the VVDS and ex- 
cludes galaxies relatively blue in R-I (i? — / < 0.7) but still 
red in B-R {B — R> 0.6), resulting in a larger correlation 
length. 



5. Discussion: evolution of the clustering length 
ro(z) from z ~ 2 

The evolution of the clustering length ro(z) from the 
VVDS first epoch "Deep" sample of galaxies with 17.5 < 
Iab < 24 is presented in Figure [7| The redshift of each 
bin is computed as the mean of the redshifts of galaxies 
in each bin. We find that ro is roughly constant or pos- 
sibly slightly increasing, within our measurement errors, 
as redshift increases, with a low value of ro — 2.35^q'^y 
h-i Mpc for z = [0.2,0.5], to ro = 3.03toje h'^ Mpc for 
z [1.3,2.1]. 

Ideally, one would like to follow the evolution of the 
clustering of the mass in the universe, translated into 
the gravitational growth of structures. To access to this 
measurement from the correlation length of galaxies re- 
quires understanding of the evolution of galaxies as com- 
plex physical systems, and to be able to relate the pop- 
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Fig. 3. The correlation function from the VVDS-02h data up to z ~ 2: (left) Wp(rp) measured in 6 redshift bins, the 
number of galaxies in each bin is indicated on the top right of each panel, (right) the correlation length ro{z) and 
correlation function slope 7 measured from fitting Wp{rp) assuming a power law form ^(r) = {r/ro)^'^ . 68%, 90% and 
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ulations of galaxies observed at different redshifts as de- 
scendents of a well-identified original population of matter 
halos at early epochs. Although the VVDS I-band mag- 
nitude selection is the minimum selection bias one can 
impose selecting a distant sample, this faint magnitude 
limit leads to a broad redshift coverage, and to a popula- 
tion mix in the VVDS sample that changes with redshift 
as described in Section 2.2. Over the redshift range con- 
sidered here, the luminosity function of galaxies in the 



VVDS is strongly evolving up to ~ 1.5 — 2 magnitudes 
llllbert et alj.l2005l) . The interpretation of the evolution of 
ro(z) is therefore not direct. 

In the lower redshift bin z < 0.5, we find a correla- 
tion length smaller than any other population of galaxies 
observed today. The mean absolute luminosity of the low 
redshift sample is Mbab = —17.5 with a significant num- 
ber of objects fainter than Mb^b — ~16 (Figure|21), while 
bright galaxies are under-represented due to the small vol- 
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measured in the DEEP2 to ease rp comparison. 

ume available at these low redshifts. Not surprisingly given 
the faint Tab — 24 cutoff in the VVDS, this makes it the 
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Table 1. Measurements of the correlation length ro(z) and the correlation function slope 7(2). The associated la 
errors are reported, vq values are computed both with letting 7 free, and setting 7 to the m ean 1.78 in the r ange 
z=[0.2,1.3] in the VVDS-F02 and VVDS-CDFS. To compare to the DEEP2 measurements of ICoil et all l|2004 . we 
have set 7 = 1.66 for the VVDS-F02-DEEP2 in z=[0. 7,0.9] and z=[0.9,1.35], and 7 = 1.68 for the VVDS-F02-DEEP2 
in z=[0.7,1.35]. The redshift range, number of galaxies used, mean absolute magnitude Mbab ^'^d I'cst frame color 
B — Iab{z — 0), are indicated for each subsample. 
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faintest galaxy population for which 3D clustering has ever 
been probed at low redshifts, about 1.5 magnitude fainter 
than the 2dFGRS sample wi th a mean Mh = - 17.98 
equivalent to Mb^b = -18.9 l|Norberg et ali l200ll) . This 
low clustering of the low luminosity population of galaxies 
measured in the VVDS is roughly consistent with an ex- 
trapolation to fainter luminosities of the trend for a lower 
correlation length as galaxies become fain ter, as observed 
locally in the 2dFGRS an d SDSS surveys ijNorberg et all 
l200ltlZehavi et a l.'.'2002^. and at intermediate redshifts in 
the CN0C2 survey (ShcDhcrd ct al., 2001). 

As redshift increases, the VVDS probes more of the 
bulk of the general population of galaxies. Galaxies sam- 
pled at redshifts 0.5 — 1.1 become increasingly similar in 
luminosity and color to the population sampled by the low 
redshift 2dFGRS and SDSS surveys. As can be seen from 
Figure [3 at z ~ 1 we measure the clustering of galaxies 
with a mean absolute magnitude Mb^b — ~18.5, after 
taking into ac count 1 m agnitude of luminosity bright- 
ening at z ~ 1 ijllbert et al. . 2005), thus comparable to the 
bulk of the galaxies probed by 2dF and SDSS at z ~ 0.1. 

Over the redshift range z=[0.5,l.l], we find a slightly 
increasing correlation length tq ~ 2.2 — 2.9 h~^ Mpc. For a 
similar, blue-selected, population of galax ies at z cii 0.15 , 
the 2 dFGRS finds vq ~ 4.3-4.6 h'^ Mpc l)Norberg et all 
l200l|) . using the same technique we use here. Note that 
at a redshift similar to the 2dFGRS, the SDSS tends to 
sample a different mix of morphological types, due to its 
red-based selection, and for this reason measures a larger 



correlation length, rp = 6.1 4±0.18 h^^ Mpc ijZehavi et alJ 
12002); see ..Haw k ins et al.l |20Q3) for discussion). Our re- 
sults therefore seem to indicate that the amplitude of the 
correlation function of galaxies which would have lumi- 
nosities AIbab — —19-5 after correction for luminosity 
evolution is about 2.5 times lower at z = 1 than observed 
locally by the 2dFGRS. 

While we know that galaxies with similar luminosities 
at different redshifts may not b e tracing the underlyin g 
dark matter in a similar way fsee lMarinoni et all l)2005|) ). 
this result is in qualitative agreement with the expecta- 
tions of the simple gravitatio nal growth of primordia l fluc- 
tuations (see e.g. Figure 5 of I Weinberg et all l)2004|) '). 

In the highest redshift bins z=[l.l,1.3] and z=[1.3,2.1], 
we are measuring the correlation function of the brightest 
Mb^b ^ —19.5 galaxies, with a mean of Mb^b ~ —20.5 
and ^ —21.0 respectively. We observe that the correla- 
tion length increases compared to the measurements in 
the range z=[0.5,l.l], up to rp ~ 3 h~^ Mpc at z=1.2 and 
ro 3.6 h~^ Mpc in the highest redshift bin to z ~ 2. The 
luminosity function of galaxies at these redshifts shows a 
marked evolution, equivalent to an increase in luminosity 
of ~ 1.5 and ^ 2 magnitudes in z= [1.1, 1.3] and z=:[1.3,2.1], 
respectively l|llbert et alll2005|) . These galaxies will there- 
fore be expected to have a mean Mb^b ~ —19.5 at low 
redshifts after evolution. Again, comparing the clustering 
observed in the VVDS to the clustering length of ro ~ 5 
h-^ Mpc measured for Mb^b = -19-5 in the 2dFGRS 
l|Norberg et al.„ 200l| . we find that the clustering ampli- 
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tude has increased by a factor ~ 2.4 from z ~ 1.3 — 1.5 to 
z = 0. 

The evolution of the correlation length observed in our 
data is in broad agreement with the results of computer 
simulations of galaxy formation and evolution. In their 
SPH simulation. IWeinberg et all l)2004 find that the clus- 
tering length of galaxies decreases from tq ~ 4.2 h~^ Mpc 
at z = 3 to a minimum ro 3.0 h~^ Mpc at z — 1.5, 
then increases again to ro ~ 4.0 h~^ Mpc at z = 1. 
Note, however, that these predictions refer to the correla- 
tion length of the same class of galaxies, corresponding to 
M > 5 X 1O^°M0, ideally followed at different redshifts. 
In a real, magnitude-limited observation, as we have dis- 
cussed, this effect is related to the different range of lumi- 
nosities sampled at each redshift, and the changing clus- 
tering strength at different luminosities. Besides the sim- 
ulated clustering of galaxies, the dark matter correlation 
length is expected from theory and N-body simulations 
to drop steeply with increasing redshift, from ro '--^ 5 h~^ 
Mpc at z=0 to rn ^ 1.8 h~^ Mpc sX z — 1.5 (see e.g. 
IWeinberg et"!!! l|2004j) l. In comparison to these predic- 
tions, our data clearly show that the clustering evolution 
of galaxies does not follow the predicted trend for dark 
matter. 

6. Conclusions 

We have computed the evolution of the correlation func- 
tion ^(rp, tt) and its integral along the line of sight Wp{rp), 
from the VVDS first epoch "deep" survey. The VVDS 
contains a large spectroscopically-selected sample of 7155 
galaxies representative of the global galaxy population in 
the redshift range z=[0,2.1], in a large O.Gldeg^ total area. 
The correlation length tq is observed to be low, vq ~ 2.2 
h~^ Mpc, for the low redshift z < 0.5 population, indi- 
cating the low clustering of the very low luminosity pop- 
ulation sampled in this redshift range. Over the redshift 
range z=[0.5,l.l], the correlation length of the population 
of galaxies, with a luminosity range comparable to the 
lower redshift 2dFGRS and SDSS, stays roughly constant 
with tq ~ 2.8 h~^ Mpc. At the highest redshifts probed in 
this paper, z=[l. 1,2.1], we find that the correlation length 
increases slightly to tq ~ 3.6 h"-'^ Mpc. 

After applying the same selection function as in the 
DEEP2 survey, our results are found to be i n excellent 
agreement with the results of ICoil et alJ l(2004l) . However, 
the significantly different DEEP2 selection function ex- 
cludes up to 41% of the galaxy population observed by the 
VVDS, making the correlation length of the DEEP2 sam- 
ple larger than in the purely magnitude-selected VVDS. 
The VVDS results demonstrate that the clustering of the 
whole galaxy population at the same redshift is indeed sig- 
nificantly lower, with tq ~ 2.9 in the VVDS vs. tq ~ 3.5 
h~^ Mpc in the DEEP2, further emphasizing the impor- 
tance of our simple, purely magnitude-limited selection 
function. 

Our measurements clearly show that the correlation 
length evolves only slowly with redshift in the range 



0.5 < z < 2, in a magnitude limited sample with 17.5 < 
Iab < 24. Taking into account the different VVDS galaxy 
populations probed as a function of redshift, with intrin- 
sically brighter galaxies probed as redshift increases, we 
find that the clustering of galaxies at z^ 1 — 2 in the 
VVDS is about 2.5 times lower in amplitude than for the 
galaxies probed by the 2dFGRS at z~ 0.15, for popula- 
tions with similar absolute Mb magnitudes. This result 
provides unambiguous evidence for clustering evolution. 

Our results are in broad agreement with simula- 
tions account ing for both gravitational growth arid bary - 
onic physics l|Weinberg et all l2004l [Benson et all kOoA . 
These simulations show that the underlying dark matter 
correlation evolves strongly with redshift, as expected in 
a hierarchical growth of structures. Our observation that 
the clustering of galaxies does not follow such a strong 
evolution therefore fully supports the model in which lu- 
minous galaxies at z = 1 — 2 (or 9 — 10 billion years ago) 
trace the emerging peaks of the large-scale dark-matter 
distribution and implies a strongly evolving galaxy bias. 

We will investigate the evolution of the g alaxy - dark 
matter bias elsewhere l|Marinoni et al.L l2005(l . A more de- 
tailed study of the dependance of clustering on luminosity 
and galaxy type will be presented in forthcoming papers. 
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